INT-PUB 03-15, MIT-CTP 3402 



A lattice theory for low energy fermions at nonzero chemical potential 
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We construct a lattice theory describing a system of interacting nonrelativistic spin s = i fermions 
at nonzero chemical potential. The theory is applicable whenever the interparticle separation is 
large compared to the range of the two-body potential, and does not suffer from a sign problem. 
In particular, the theory could be useful in studying the thermodynamic hmit of fermion systems 
for which the scattering length is much larger than the interparticle spacing, with applications to 
realistic atomic systems and dilute neutron gases. 
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INTRODUCTION 



In this Letter we consider dilute fermion systems with 
attractive interactions, for which the effective range for 
two-body scattering is much less than the inter-particle 
spacing. Such systems are characterized by a dimension- 
less number rj — (an^^^), where a is the two-body scat- 
tering length, and n is the density. Such systems are well 
known to exhibit fascinating nonperturbative behavior. 
For 77 small and negative (weak attraction) one finds the 
BCS solution with pairing and superconductivity. With 
77 small and positive, corresponding to strong attraction 
with a two-body bound state well below threshold, the 
bound pairs will Bose- Einstein condense. In each of these 
cases the behavior is nonperturbative, but since the ef- 
fective interaction is weak, the system can be successfully 
described in a mean field approximation. On the other 
hand, dilute fermion systems for which the parameter \r]\ 
is large are not amenable to analytical treatment. Phys- 
ical realizations include both dilute neutron gases (the 
neutron-neutron scattering length is more than an order 
of magnitude greater than its effective range) as well as 
cold, dilute gases of fermionic atoms tuned to be near 
a Feshbach resonance 0, 0, 13 ■ Recently there has been 
intense interest in exploring such systems experimentally 
Q. In the hmit that \r]\ ^ 00 one expects to see uni- 
versal behavior, so that the same dimensionless physical 
constants will apply equally to the atomic and nuclear 
systems. 

Because of the prospects of exploring ultracold 
fermionic atoms at large I77I experimentally, it is of great 
interest to understand such systems theoretically, which 
necessitates numerical studies. A recent numerical study 
of dilute fermions at large \ri\ is found in Ref. 5]. The 
authors of that work used the fixed-node Diffusion Monte 
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Carlo method to extract the parameter ^, defined as the 
energy per particle relative to the value for a nonintcr- 
acting Fermi gas. Since this approach is variational, the 
result ^ = 0.44 ± .01 obtained is an upper bound on the 
true quantity. 

The calculation 0| is performed at fixed particle num- 
ber A^, up to = 38. It is desirable to have a numer- 
ical approach which is not variational, and which can 
probe thermodynamic properties of dilute fermion gases, 
such as the critical temperature for the pairing transition. 
In this Letter we propose such a method, in which the 
fermions live on a spacetime lattice at nonzero chemical 
potential. We exploit the fact that all relevant scatter- 
ing information at low density is contained in the scat- 
tering length, so that the complicated two-body inter- 
actions can all be replaced by an effective field theory 
with a single attractive four-fermion interaction. While 
lattice formulations of interacting fermions at finite den- 
sity typically encounter a "sign problem" , we show that 
is not the case here. There is another issue that requires 
attention, however, relating to possible instabilities in- 
troduced by a purely attractive interaction. We perform 
an analytical calculation that suggests such instabilities 
may be avoided. 



II. THE LATTICE FORMULATION OF THE 
EFFECTIVE THEORY 

To construct a lattice version of the continuum the- 
ory, we first analytically continue to Euclidean spacetime, 
represented by an A^,- x A^^ lattice, where A^,- and Ng are 
the number of sites in the time and space directions re- 
spectively. We set h — 1 and measure all dimensionful 
quantities in units of the lattice spacing, which is taken 
to be the same in the space and time directions. We 
impose periodic boundary conditions in the space direc- 
tions and anti-periodic boundary conditions in the time 
direction. This allows the Euclidean path integral to be 
interpreted as the finite temperature partition function, 
with inverse temperature P = N-r- The fermion fields at 
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site n are denoted by f/'n and 'i/'n, and are independent, 
two-component Grassmann spinors. 

The continuum Euclidean action for free 
fermions at nonzero chemical potential is 
/ d'^x ip [dr ~ V^/2AI — ^) -ifj. The corresponding 

lattice action is Sq = V'n(-K^oV')n, where 

(ifoV'jn = [ipn - e'^V'n-eo) - 2^ ^^^^^ 

1=1 

. (1) 

In these expressions is a unit vector in the a direction 
(a = corresponding to Euclidean time r), M is the 
fermion mass, and /i is the chemical potential which is 
treated like the time component of an imaginary gauge 
field. This follows the prescription in ref. (except that 
we have the opposite sign for fi) and ensures that the free 
energy does not have spurious cutoff dependence. 
The free propagator computed from 5*0 is 



Go(r,p) = 



1 if r > 



where 



(2) 



1 + Ap/M 



A 



p^2^sin^|. (3) 



The time t is an integer, corresponding to propagation 
by r lattice spacings in the time direction, while the mo- 
mentum components are given by — Pij^, with pi 
being integers in the range ~Ns/2 < pi < Ns/2. 

Note that for < 0, ^p < 1 and in the limit — > oo 
(zero temperature) — > and there is no propagation 
backward in time, 



Go(t,p) 



fr 



M<o,w,^oo (1 + Ap/M) 



(4) 



indicating the absence of antiparticles or holes in this 
nonrelativistic theory at zero density. 

Fermion interactions can be represented by an effective 
field theory with a four-fermion interaction the strength 
of which is tuned to reproduce the physical two-body 
scattering length. On our lattice we generate the inter- 
action by means of a non-propagating scalar auxiliary 
field if coupled to the fermions. By situating ip along 
time-links of the lattice we can eliminate fermion loops 
lying on surfaces of constant Euclidean time, simplifying 
the analysis. Since one is interested in pairing correla- 
tions, it is convenient to also introduce a constant com- 
plex source J for fermion pairs so that one can study the 
pairing transition at finite volume. These considerations 
lead us to the action 

S = Y1 [V'n(ifV')n + ^"iVn + \ iJi^nCr2^n + h-C.) 
n 

(5) 



An important property of this theory is that after in- 
tegration over ip and tp, the remaining integral over (p 
has positive semidefinite measure, so that Monte Carlo 
integration methods may be applied. This result is easy 
to see if the source J is neglected, in which case inte- 
grating out the fermions yields a factor of det K. K is 
trivial in spin-space, and can be written as K = K Is 
where l^. is the two dimensional identity matrix acting 
on spinor indices, and K acts only in coordinate space. 
Thus detK = {det K)^ which is positive semidefinite, 
since if is a real operator. 

Including the constant source J, the fermionic action 
may be rewritten as ^'i'^A'if , where we have defined 



* = 



T > 



A = 



-iJ A'T 

K ~i.r 



i<72 
10-2 



(6) 



Each block in this matrix is 2N dimensional, where N 
is the number of lattice sites; K acts trivially on spinor 
indices, the Pauli matrix cr2 acts trivially on coordinate 
indices, and J is proportional to the identity matrix in 
both spaces. Note that — K'^ since K is real. Inte- 
grating out the fermions yields a purely bosonic theory, 
where the path integral over the Lp field is weighted by 
exp{—^ vii) times the pfaffian of A: 



Vi[A] 



-iJ ift 

k -i.r 



<^ \J\' 



(7) 



Each block in the above expression is A^-dimensional. 
Thus Pf[-4] is real and positive semidefinite for all real 
values of J, /i, and (/3„, up to the irrelevant constant sign 
a = (—1)^, and there is no sign problem encountered 
with this lattice formulation. Similar analyses exist for 
two flavor QCD and the NJL model [3. 

We have introduced in Eq. ^ a new parameter m^, 
which determines the strength of the two-body interac- 
tion. This parameter can be related to the two-body 
scattering length a by summing the ladder diagrams with 
Lp exchange between the fermions at zero external mo- 
menta, zero temperature, infinite volume, and zero chem- 
ical potential, /i ^ 0~; the result may be equated to 
(see (^]). The result is 



m 



1 

Ana 



L{M) , 



(8) 



where L{M) is given by the integral 



L{M) 



1 f (fp 



(9) 



2 J (27r)3 Ap + A2 /2M 



where the integral is over the Brillouin zone, \pi\ < tt. 

Note that since must be positive, scattering lengths 
satisfying i > AttL^M) are inaccessible in our theory 
without introducing an imaginary coupling for (p and sac- 
rificing positivity of the measure. These values for the 
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scattering length correspond to two-body bound states 
with 0{1/M) binding energy in lattice units, which are 
not of physical interest in any case. 



III. PHASE STRUCTURE OF THE LATTICE 
THEORY 



An interacting continuum limit of our lattice theory 
requires that we be able to find a point in the phase dia- 
gram where both the interparticle spacing n~^/^ and the 
scattering length a diverge in lattice units. It does not 
require that the mass M vanish, however. From Eq. © 
we see it is always possible to tune the parameter m? so 
that a diverges. What is less obvious is that we can then 
choose /i so that the particle density n becomes arbi- 
trarily small. In general one would expect to find a line 
of first order phase transitions in the (Af, i) plane for 
sufficiently large values of ^ , corresponding to a strongly 
attractive four-fermion interaction. This follows from the 
well-known result in the continuum that a fermion sys- 
tem at fixed chemical potential with a purely attractive 
interaction has no ground state, being unstable against 
the formation of infinitely dense matter. On the lattice, 
where the maximum attainable density is limited, one 
would expect that at = the ground state of the sys- 
tem would jump from n = to n = 0(1) in lattice units 
for a sufficiently strong attraction. The existence of a 
continuum limit then becomes whether for any mass M 
the n — Q phase at /i = persists when 1/a — > 0. 

Therefore the first thing to compute in lattice theory 
is the relation between the chemical potential /i and the 
fermion density n. For n we take 



(10) 



where 'P(^) (^e""^ Pf [yl]^ , and Pf[^] is given 

in Eq. 10. This definition of n is more convenient than 
din Z / djj.^ and the two definitions agree when either ex- 
pression is small. 

To compute the density n requires a full-scale un- 
quenched simulation, which is beyond the scope of this 
paper. However it is instructive to determine n in the 
leading semi-classical (mean-field) approximation. We 
compute the effective potential obtained by integrating 
out the fermions, assuming a classical background {ip) 
field; we then determine {(p) to be the value that mini- 
mizes the effective potential. 

Setting J = in Eq. iQ (and dropping the sign <j) we 
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FIG. 1: The density phase diagram in the (M, a~^) plane 
at /I = in the mean-field approximation for M > 0.1. The 
curve corresponds to a first order phase transition in the ex- 
pectation value of the fermion number density n, with n = 
in the lower region and n — 0(1) in the upper region. The de- 
sired continuum physics only exists in the lower phase. The 
black region in the upper left corner is inaccessible for real 
scalar-fermion coupling. 



can compute Pf[.A] exactly for constant ip: 

Ff[A] = detK^'K 

= n [(1 + Ap/A/)^- + e^^^l + pf^ 
{pT 



E 



In 1 



M 



+ 



ln(l + (Cp(l + y))^-) " 



(11) 

In the Nt-,Ns — > oo limit (large volume, zero tempera- 
ture) this becomes 



Vi(^) 



BZ 



(2^ 



ln(l + Ap/A/) 



+ ln(Cp(l + ¥'))e(ep(l + ^)-l) 



The momentum integration in the above expression is 
over the Brillouin zone, pi £ (— tt, tt]. The first term in 
Vi is just the contribution from free fermions; the second 
term is the quantum correction to the tree level potential, 
Vo{(p) = ^m?'pP'. In the mean-field approximation the 
expectation value {(p) is given by the minimum of the 
effective potential Veff(<p) = (Vq + Vi). 

The effective potential yields the ^ — Q phase diagram 
displayed in Fig. ^ The horizontal axis is the fermion 
mass Af, while the vertical axis is the inverse scattering 
length, . Moving up the vertical axis corresponds 
to increasing the strength of the two-body attraction, 
and > describes systems with a two-body bound 
state. The curve corresponds to a jump in the fermion 
density from zero in the lower region, to 0(1) in the 
upper region. We see that there is a region far from 
the phase transition for M ^ 1 that allows small of 
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FIG. 2: Density n as a function of ^ for M = 1 and a oo. 
This figure exhibits a first order phase transition as a function 
of /i. The small fi phase is where a continuum limit is defined. 

either sign. Thus one may expect to be able to perform 
simulations of fermions at finite density, where both the 
interparticle spacing n~^/^ and scattering length \a\ are 
large compared to the lattice spacing, without requiring 
1 77 1 to be small. 

If one were to plot the chemical potential fi along a 
third axis of the phase diagram, then the critical line 
we have displayed would become a critical surface. In 
Fig. 121 we show the fermion number density n plotted as 
a function of /x, for fixed M = 1 and = 0. One 
sees that there is a regime where n grows slowly with 
/i, but that at a critical chemical potential, the fermion 
density jumps discontinuously as the trajectory crosses 
the critical surface. We see that at least in the mean- 
field approximation, there is a large region in parameter 
space where we can define continuum theories with any 
value oi rj = (an^/^) we desire. 

One of the first tasks of a fuU-fiedged simulation will be 
to determine the analog of Fig. Q] beyond the mean field 
approximation by evaluating the expression in Eq. 11U|) 
and to show that a region in parameter space with an 
interesting continuum limit persists. An interesting cal- 
culation to pursue subsequently would be to map out 



the critical temperature for the superfluidity phase tran- 
sition with order parameter a function of the 
dimensionless quantity rj. To do this one would need to 
simulate the system at finite source J, and extrapolate 
to the J = 0, infinite volume limit. 



IV. DISCUSSION 

There are numerous ways in which the lattice theory 
we have presented may be modified or extended without 
sacrificing its crucial property of positivity. Several ob- 
vious possibilities are: (i) to include a kinetic term for 
the field, allowing for attractive fermion interactions 
extended in space (for example, a Yukawa interactions 
for the simulation of neutron matter, for which p-wave 
interactions are expected to be important ^3); (ii) 
</5 field could be given derivative and spin-dependent cou- 
plings, perhaps allowing one to include the contributions 
from the two-body effective range or scattering in higher 
partial waves; (iii) more flavors can be introduced, which 
could provide an opportunity to study the importance 
of three-body contact interactions and the possibility of 
renormalization group limit cycles 0, 0, 0, 0| . 
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